Measurement of foliar H2O2 concentration can be an indicator of riparian vegetation management

Riparian vegetation is frequently exposed to abiotic stress, which generates reactive oxygen species (ROS) caused by strong differences in a river’s hydrological conditions. Among different ROS, hydrogen peroxide (H2O2) is relatively steady and can be measured appropriately. Thus, the quantification of plant H2O2 can be used as a stress indicator for riparian vegetation management. The current study examines the spatial distribution of plants by riparian vegetation communities across the elevation gradient of riparian zones through quantification of environmental stress using foliar H2O2 concentration. The trees Salix spp., Robinia pseudoacacia, Ailanthus altissima with Juglans mandshurica, and the herbs Phragmites australis, Phragmites japonica, and Miscanthus sacchariflorus were selected for this study. Leaf tissues were collected to analyze H2O2 concentration, meanwhile riparian soil was sampled to measure total nitrogen (TN), total phosphorus (TP), and moisture content. The H2O2 concentration of tree species increased with higher soil moisture content, which was negatively correlated for Salix and herb spp., in which H2O2 concentration always decreased with high soil moisture. In this study, we found a unique significant interaction between soil moisture content and H2O2 concentration, both positively or negatively correlated relationships, when compared with other parameters, such as TN or TP concentrations or TN: TP in riparian soil. The species-specific distribution zones can be explained by the H2O2 concentration in the plant for gravelly and sandy channels on a theoretical range of soil moisture. Each species’ H2O2 concentration was estimated through derived equations and is directly related to an elevation above the channel. The comparison with the observed distribution of plant elevations in the field indicated that all species showed a spatial distribution that acts as species-specific elevations where H2O2 concentrations stayed below 40 μmol/gFW. Hence, the present study suggests that foliar H2O2 concentration can be a useful benchmark for the distribution potentiality of riparian vegetation.


Results
Edaphic condition of the riparian zone. The soil moisture content gradually decreased with increased elevation away from the channel at both the Arakawa River and Hii River sites (r = −0.90, p < 0.01 for the Arakawa River and r = −0.98, p < 0.01 for the Hii River) (Fig. 1). The distributions were empirically given as follows: For the Arakawa River (gravelly) site (r = 0.99): For the Hii River (sandy) site (r = 0.99): where elevation is given in meters above the ordinary water level.   55,56 for comparison. The Salix species localized generally at higher soil moisture sites than other species [57][58][59][60][61] . Regardless of rivers and the reservoir site, the foliar H 2 O 2 concentration of all Salix species had a significant decreasing trend with increasing soil moisture content (r = −0.89, p < 0.01 for S. pierotii; r = −0.92, p < 0.1 for S. gilgiana; r = −0.5, p < 0.1 for S. subfragilis). This continued until the soil moisture content reached 35%, keeping the nearly constant value with www.nature.com/scientificreports/ higher soil moisture content. Thus, the lowest leaf H 2 O 2 concentration (20 µmol/gFW) was recorded at higher than 35% soil moisture content. Other tree species exhibited positive correlations between soil moisture content and H 2 O 2 concentration (r = 0.66, p < 0.01). The value of H 2 O 2 concentration (approximately 15 µmol/gFW) was lowest at around 5% soil moisture content (the driest condition), then rose with increasing soil moisture until approximately 20-30%, where it attained about 40 µmol/gFW. Foliar H 2 O 2 concentration of herbaceous species concerning soil moisture content is presented in Fig. 4b. Regardless of sites, similar species-specific trends were observed. M. sacchariflorus had a higher H 2 O 2 concentration at low soil moisture conditions, which decreased with increasing soil moisture until 15% (r = −0.94, p < 0.001). It slightly increased again with higher moisture content. Phragmites spp. showed a uniquely decreasing trend with soil moisture content (r = −0.84, p < 0.001). Overall, there was a negative correlation between herb species and soil moisture content (r = −0.80, p < 10 −5 ).

Foliar H 2 O 2 concentration with respect to other factors. TN and TP contents of plant biomass have
organ-specific differences. Plant's N: P ratios < 10 and > 20 often (not always) correspond to N-and P-limited biomass production in short term periods respectively, although it can vary in the long term 62,63 . The TN contents of plant biomass were approximately 0.9%, 0.4%, and 2.1% of roots, stems, and leaves of Salix spp., respectively, and 1.3%, 0.6%, 3.0% for R. pseudoacacia. TP contents were 0.1%, 0.01%, 0.12% for Salix spp. and 0.05%, 0.12%, and 0.07% for R. pseudoacacia. For herbs, whole-plant values were 1.9 ± 0.25% for TN and 0.19 ± 0.03% for TP. The nitrogen and phosphorus ratio of plant biomass's range was 10 to 30. This was much larger than the nitrogen and phosphorus ratio of soils, which was 1 to 4. Therefore, compared with phosphorus, nitrogen seems to be more restrictive for plant growth. However, given the correlation between biomass (of P. australis or M. sacchariflorus) and riparian soil, the TN and TP contents were 0.24 ± 0.05 (p = 0.37) for TN and 0.1 ± 0.37 (p = 0.71) for TP (data not shown). There was no correlation between biomass and TN or TP contents in the riparian soil.
The H 2 O 2 concentration with respect to riparian soil TN or TP contents is shown in Fig. 5      The non-linear regression analysis was performed for each species to find out effect of different parameters (soil moisture, H 2 O 2 , TN, TP, TN: TP). Salix spp., R. pseudoacacia, A. altissima, and J.mandshurica individually show significant correlation between soil moisture content and H 2 O 2 (p < 0.001, p < 0.001, p = 0.002, and p < 0.05 respectively), whereas no significant results were observed (p = 0.459 for Salix species, p = 0.884 for R. pseudoacacia, p = 0.186 for A.altissima, and p = 0.652 for J.mandshurica) among parameters (soil moisture, TN, TP, and TN:TP). Herb species also exhibit a similar type of significant trend in the nonlinear regression analysis. Phragmites spp. and M. sacchariflorus indicate significant correlation between soil moisture content and H 2 O 2 (p < 0.001 for both species) On the contrary, no significant results were noticed (p = 0.076 for Phragmites spp., and p = 0.138 for M. sacchariflorus) among parameters (soil moisture, TN, TP, and TN:TP).
We cannot find significant trends between H 2 O 2 concentration and riparian soil TN or TP. One of the reasons is because the variation range of TN and TP was too small to affect H 2 O 2 concentration in the observed sites 64  Modeling foliar H 2 O 2 concentration. From the above discussion we can conclude that due to the insignificant impact of nutrients especially TN or TP concentration, the H 2 O 2 concentration in plants is solely related to soil moisture content. The simple empirical equations can be formulated based on the Monod equations. The equations show the fundamental structure as easily explicit and broaden the underlying mechanism of the species.
For tree species: Salix spp.
(for the range of soil moisture < 45%, r = 0.92, p < 10 −9 ). Miscanthus sacchariflorus: (for the range of soil moisture < 30%, r = 0.89, p < 10 −3 ); where soil moisture is given by %. Soil moisture distribution depends on the sediment particle size composition. This relationship was similar for either gravelly or sandy rivers (Fig. 1). For simplicity, soil moisture at the ground surface at the Arakawa River site, Eq. (1), and the Hii River site, Eq. (2), are used as representative of those gravelly rivers and sandy rivers in general. Substituting Eqs. (1) and (2)  Both Salix spp. and other tree species are restricted to areas with stress levels less than 40 to 45 μmol/gFW. Therefore, Salix spp. can only distribute starting from the ordinary water level up to 4 m in gravelly rivers and 6 m in sandy rivers. In contrast, other tree species can distribute from 1 m upwards in gravelly rivers and 3 m upwards in sandy rivers. For herbaceous species, all species areas appear where soil moisture results in an H 2 O 2 concentration below 40 μmol/gFW. Phragmites spp. can, therefore, distribute at elevations up to 3 m in gravelly rivers and more than 5 m in sandy rivers. M. sacchariflorus grows at higher elevations and cannot grow below1.5 m above the ordinary water level.
H 2 O 2 is generated under stressed conditions. Thus, the intensity of stresses is correlated to the H 2 O 2 concentration. H 2 O 2 concentrations of light-exposed and dark-adapted samples are shown in Fig. 8. There was no difference in H 2 O 2 concentration between the two treated leaves (p = 0.96) for all tested species, indicating that dark conditions do not cause significant stress.
To ensure the efficacy of our derived equation, we try to analyze their suitability with the general additive model (GAM). We found that GAM supports our derived equations appropriately among different parameters

Discussion
In Japanese rivers, the volume of sediments has substantially decreased over the last five decades because of the construction of dams upstream, gravel mining, and the reduction of sediment inflow due to the afforestation of the upstream mountainous areas. Therefore, the riparian morphology and sediment characteristics remain relatively stable, even though it is subjected to recurrent large flooding events 17,27,65,66 . Most of the sampling points were inundated only during a large flooding event, which occurs approximately once every 10 years, with the inundation lasting for a maximum of one week. Therefore, the measured soil moisture content is characteristic of the study site most of the time. The moisture content that the riparian zone reaches is highly affected by sediment particle composition; it is higher at sandy channels than at gravelly channels, as shown in the present study. There were good agreements between the results of previously observed cases (from Eqs. (1) and (2)) and the current study at both river sites. Several other factors can affect riparian vegetation stress. In the riparian zone of the observed reach, generally nitrogen rather than phosphorus becomes critical for plant growth 17,67-69 , and it may thus influence the growth of vegetation via soil moisture. However, in the present study, there was no correlation between plant tissue H 2 O 2 concentration and TN along with TP concentration in the riparian soil for the tree and herb species as shown in Fig. 6. Therefore, nitrogen and phosphorus were considered as neither major factors for limiting growth nor regulating species distribution in both rivers. Also, solar radiation did not affect the H 2 O 2 concentration as expressed in Fig. 5. The last major flooding was more than two years ago, and no erosion or deposition trace was observed on the soil surface 26 . Thus, mechanical disturbances on the ground surface during floods likely did not affect the distribution of plants. Therefore, soil moisture seems to be the major component of stress for the plants in this case.
Threshold H 2 O 2 concentration . When a plant is exposed to environmental stress, it generates H 2 O 2 in different organelles. This occurs partly in the scavenging process of various reactive oxygen species and partly directly 37,38 . However, some H 2 O 2 is scavenged by the antioxidant activities and finally converted to water and oxygen. Under high-stress conditions, plant tissues become damaged and die off due to the high concentration of H 2 O 2 . Therefore, there is likely a threshold concentration below which plants do not deteriorate or die 47,48 . In the present study, no cases with H 2 O 2 concentrations higher than 40-50 µmol/gFW were observed for tree species. Concentrations above 40 to 50 µmol/gFW of H 2 O 2 seem to be too high to maintain a healthy condition for a tree species. For example, Salix subfragilis in the Miharu Dam highly deteriorated at H 2 O 2 concentration levels ~ 40 to 50 µmol/gFW, which was caused by at least 5 months of low soil moisture 44 . The low soil moisture condition seems to be a high stressor for Salix spp. At the same time, each species distributed only the elevation where the H 2 O 2 concentration becomes lower than these values.
In the present study, R. pseudoacacia and A. altissima growing on low elevation sites along the Arakawa River had a high H 2 O 2 concentration of 40 to 50 µmol/gFW. These plants may experience greater stress with high soil moisture. All R. pseudoacacia and A. altissima plants intrude and grow at the lower elevation sites (below 1 m, 30% moisture; Fig. 1), had died off over the course of several years 26,52,54 without being subjected to flood disturbance. This situation indicates that more than 40-50 μmol/gFW of H 2 O 2 concentration reflects harsh conditions for R. pseudoacacia and A. altissima (Fig. 2a). In contrast, all S. gilgiana growing at elevations above 3 m had disappeared after several years 26,52,54 . S. gilgiana cannot spread at this elevation, probably because of very low soil moisture content. For herbaceous species, no Phragmites spp., were found at elevations with less than 10% soil moisture in the gravelly channel and 20% in the sandy channel. This corresponds to an H 2 O 2 concentration of 35 μmol/gFW and above. M. sacchariflorus distributed at relatively low soil moisture, with the lowest being about 5%, corresponding to higher than 40 μmol/gFW of H 2 O 2 concentration. Although the water potential was not measured in this study, the sediment's logarithmic scale of suction (pF) was approximately 4 at 10% moisture content 69,70 . The pF ~ 3 to 4 corresponds to the wilting point and is nearly the critical condition for the continued growth of many herbaceous plants 71,72 . The 40-50 μmol/gFW of H 2 O 2 concentration obtained for the existing plants seems to be the critical value for survival. The lowest elevation where M. sacchariflorus was found was approximately 1.5 m above the ordinary water level, where the plant had about 30 μmol/gFW of H 2 O 2 concentration. The location seemed close to the lowest elevation for the distribution, although this was not absolutely critical. Thus, the H 2 O 2 concentration of M. sacchariflorus is U-shaped with respect to soil moisture. For other species, Salix spp. and Phragmites spp. were not found in water, and the H 2 O 2 concentration is unknown, and Characteristics of tree species. Salix species usually are dominant at low elevation zones [57][58][59][60][61] . In the present study, Salix species always developed above the ordinary water level, and no Salix was found in the constantly inundated zone. Salix spp. make use of hydrochory: dispersal of seeds through water, which encroaches onto the shoreline in spring 44,58 . Therefore, the distribution of Salix spp. seems to reflect the physiological survival of occupied seeds.
For Salix spp., H 2 O 2 concentration decreases accordance with soil moisture content till 40%. This situation indicates low stress in high soil moisture content. Then, it slightly increased with more than 40% soil moisture. The preferable soil moisture of these Salix species appears to be around 40%, which was near saturation. Several experiments indicated that even steady inundation does not lead to high stress as long as oxygen concentration is sufficiently present in the substrate 61 . Other tree species in the present study provided a relatively unique increasing trend of H 2 O 2 concentration from 5% of soil moisture to higher levels. The soil moisture thus is the major component determining the potential locations for spatial distribution in the riparian zone.
For R. pseudoacacia, the highest photosynthesis rate was reported at 17% soil moisture. It was substantially lower at 8% and 24% 73 , and the highest reported leaf water potential was at 13% soil moisture, compared to lower soil moisture contents 74 . These results suggest the species' preferred specific soil moisture conditions, which is confirmed in the present results (Figs. 2a and 3a). The soil moisture content at the root zone of trees was not measured in the present study. However, it can be derived from the soil moisture distribution in Fig. 1. The soil moisture at the deepest root zone (i.e., 4 m above the ordinary water surface and 2.5 m below the ground surface) was 16%, although the surface soil moisture content was ~ 5% 52 . This species does not seem to prefer a high moisture level, although it is typically found in riparian zones. Specifically, this species can thrive at higher elevations in the riparian zone 54,75 . Since this species uses hydrochory like the Salix spp., these results also imply that seeds can reach the preferred elevations only during very high water levels. A. altissima showed a similar trend of H 2 O 2 concentration in response to soil moisture content as found in R. pseudoacacia. The H 2 O 2 concentration level of J. mandshurica was slightly lower than R. pseudoacacia and A. altissima, indicating that this species prefers higher moisture content. Its high buoyancy can be more easily dispersed downstream and captured ashore at the riverside, where the soil moisture content is more congenial to its growth.

Characteristics of herbaceous species. H 2 O 2 concentration in herbaceous species decreased with
increasing soil moisture, and this relationship was consistent regardless of species. M. sacchariflorus is one of the most common species of the riparian zone and is distributed at higher elevations than Phragmites spp. However, the H 2 O 2 concentration of M. sacchariflorus was as high as 40 μmol/gFW when it grew on soil with approximately 5% moisture content. No distribution was found at higher elevations where the riparian soil has less than the soil moisture content. Phragmites spp., on the other hand, is typically distributed in the lower elevation zone where the soil moisture content is > 10% 16 . The tissue H 2 O 2 concentration of Phragmites spp. in response to soil moisture content also indicates its preference for higher moisture conditions. H 2 O 2 concentration as the monitoring system. Riparian plants are subject to various stress factors, and they establish successfully only in areas where stress levels are sufficiently low. The riparian zone is often affected by flood disturbance, which changes the channel topography. However, although the reduction of sediment flux erodes the channel bed, it provides a stable habitat for riparian vegetation 65,66 . Groundwater level variation takes a longer period than surface water and does not change much during a short period. If the soil moisture distribution is properly maintained with an appropriate period of low flow, regardless of short flood occurrences, the elevations where these species distribute or disappear can be predicted. Measurement of H 2 O 2 can be a method to evaluate the threshold values of certain determining abiotic variables beyond which plants cannot survive. Therefore, their chance of disappearance under these conditions can be forecasted. This method is also more precise than most previous methods based on plant traits and empirical monitoring of their spatial distribution only. It is difficult to understand their preferred habitats by the general condition of their spatial distribution, owing to various types of stresses combined with different intensities and frequencies, and exposure periods 10,13,15 . Also, it is highly dependent on confounding factors like timing and unusual events. The physiological condition of each species determines their unique distribution, excluding physical mortality.

Conclusion
Riparian vegetation encounters different environmental stresses. The findings of this study demonstrate that the spatial distribution of different species in a riparian zone occurs in the context of elevation and soil moisture content. The nutrient level especially TN and TP cannot play significant role due to a small variation of TN and TP concentration. The H 2 O 2 concentration and soil moisture content solely significant in the observed sight. The species-specific distribution zones can be explained by the H 2 O 2 concentration in the plant, which could be a rapid, efficient, and reliable monitoring indicator for vegetation distribution. This study suggests that foliar H 2 O 2 concentration represents the sum of all abiotic stress, and the plant population decreases or cannot survive when it is beyond a critical value. Study site. Field sampling was conducted in the gravelly sites of the Arakawa River and the sandy sites of the Hii River to compare differences depending on sediment particle size (Fig. 9). Both sites were in the anastomosing section. The Arakawa River is located at the center of Japan, originating from the Chichibu Mountains, flowing over 173 km before draining into Tokyo Bay. The middle section of the river consists of gravel channels typical of Japanese rivers. Sampling was conducted at the Kumagaya Oaso gravel bar, 84 km upstream of the river mouth (36°8′20″ N, 139°20′35″ E). The channel slope was approximately 1/500, and the mean sediment particle size at the sampling site was approximately D 50 = 50 mm (D 25 = 10 mm, and D 75 = 99 mm). The land use in the basin upstream of the sampling site was mostly agriculture and forests. Historically, the riparian zone was occupied by non-vegetated gravel areas until 50 years ago. However, vegetation coverage has gradually increased since then due to river regulation, and now half of the area is covered with woody or herbaceous plant species. More details of the Arakawa sampling site are given in other works of the authors 16,17,26,49,52 .
The Hii River is located in the Shimane prefecture (Western region of Japan). The river originates from the Chugoku Mountains and flows first into Shinji Lake, then into Naka Lake before draining into the Sea of Japan. It extends about 153 km and drains a basin of about 2017 km 2 . Eighty percent of the basin area is covered with forest, and the downstream 20% is covered chiefly with rice fields and residential areas. The sampling site was located 14 km upstream of Shinji Lake (35°21′57″ N, 132°47′20″ E), where the channel slope was approximately 1/600, and the bed was entirely composed of sand, D 50 = 0.5 mm (D 25 = 0.1 mm and D 75 = 0.8 mm). This sand had been deposited in the past after iron mining activities upstream.
The riparian zones of both sampling sites are characterized by an upward elevation gradient across the river channel, reaching about 10 m in elevation above the water level in its ordinary flow condition. Groundwater levels have been relatively constant over the years and are tied to the ordinary river water level, at least at the locations where sampling was conducted (The Arakawa Upstream River Management Office, The Izumo River Management Office). Therefore, the riparian topsoil-to-groundwater distance increases as the distance from the river channel increases, while the overall soil moisture level decreases. The sediment's rhizospheric zone was not saturated with water, and oxygen could penetrate sufficiently except in the zone close to channels 76,77 . All sampling points were taken on the riparian slopes where the ground surface was moderately covered with vegetation without being completely shaded. Both rivers and riparian zones are subject to intense flooding events. These floods are generally short and are caused by heavy rainfall. Once every 10 years, a major flood occurs. In all cases, however, the water level returns relatively quickly to the normal condition, that is, within less than a week. Except for these events, the water surface level is relatively stable. The soil moisture pattern in the riparian zone, which decreases with elevation from the ordinary water level, has not changed significantly for decades, as evapotranspiration rate balances with precipitation rate in most seasons. These conditions make the selected field sites excellent for studying the relationship between soil moisture and plant stress.
Plant and soil sampling. Sampling was conducted between 10:00 and 14:00 on days with normal, fine weather conditions for the respective sites and more than a month after the most recent flood. In the Arakawa River, the sampling dates were May 5 and June 6, 2017, nearly 2 years after a large flood on September 9, 2015, which submerged sampling points. In the Hii River, sampling dates were October 11-13, 2016, more than 3 years after a large flood. The sampling dates were more than a year later than flooding in the respective rivers in all cases.
Dominant plant species of each site were sampled for this study. The major tree species at the Arakawa site were S. gilgiana, S. subfragilis, R. pseudoacacia, and A. altissima; the major herbaceous species were P. australis, www.nature.com/scientificreports/ P. japonica, and M. sacchariflorus. At the Hii River site, the major tree species were S. pierotii and J. mandshurica; the major herbaceous species were P. australis and M. sacchariflorus. After surveying 30 km of the riparian zone along the river reaches, three healthy mature individuals per tree species and five well-developed plants per herb species across the elevation gradient of the site were selected and marked. Selected tree species were 3 to 5 m high, and herb species were about 1.5 to 2 m high. The average root depth of the sampled plants was approximately 2 ± 0.5 m for the tree species and 0.3 ± 0.1 m for the herb species (not sampled, data based on the works) 17,26,49,50,52 .
Well-grown leaves exposed to solar radiation were carefully sampled by hand. To measure the difference between the solar radiation exposed and dark-adapted samples, additional leaves of the tree species were shaded for 30 min by covering them with a well-ventilated paper bag before sampling. All sampled leaves were sealed in a plastic bag immediately after sampling and stocked in an icebox filled with dry ice (~ − 70 °C) to transport to the laboratory, where samples were frozen at − 80 °C until analysis. It took 8 h at most to transport to the laboratory. The half-life period of H 2 O 2 is about 1.4 to 58 h in water 78,79 .
The low temperature in the stocked icebox may be a stress to changing H 2 O 2 concentration. However, the temperature declined quickly to below − 10 °C, which is lower than the freezing point of H 2 O 2 as the tissue thickness was less than 1 mm. There was no significant difference in the H 2 O 2 tissue concentrations between the icebox-stocked and immediately analyzed samples (p = 0.76) and between icebox-stocked and non-stocked samples (p = 0.90).
The exact location of all sampled plants was recorded using GPS (Garmin eTrex). Absolute elevations of sampling points and the ordinary water level above the sea level (referenced to Tokyo peil (T.P.)) were obtained from the surveyed sectional map of the Upstream Arakawa River Management Office for the Arakawa river site and the River Management Offices of Izumo for the Hii River site. Soil samples were collected in triplicate from 20 cm below the ground surface from all plant sampling points. The riparian soil was stocked in tightly sealed plastic bags and brought to the laboratory for analysis.

H 2 O 2 assays of plant leaves.
In the analytical process of H 2 O 2 concentration, the fresh plant leaves were dry frozen with liquid nitrogen and ground (~ 500 mg) together with an ice-cold 50 mM phosphate buffer (pH 6.0). Polyvinylpyrrolidone (PVP) was added to this extraction to mask the effect of phenolic compounds in the plant materials. The extraction was centrifuged at 5000 rpm at 4 °C for 15 min, and the supernatant was separated. An aliquot of 750 μL extracted supernatant was mixed with 2.5 mL of 0.1% titanium sulfate in 20% (v/v) H 2 SO 4 80 . The mixture was centrifuged at 5000 × g at 20 °C for 15 min. The intensity of the yellow color developed through the reaction was measured spectrophotometrically (UV mini 1210, Shimadzu, Japan) at a wavelength of 410 nm. The absorption at 410 nm includes the effect of other soluble compounds 81 . Thus, the H 2 O 2 concentration was calculated from the slope of the standard curve obtained from the known H 2 O 2 concentration, which was offset, derived by the intercept absorption rate with zero H 2 O 2 concentration samples 82 . The results were compared with those of the e-FOX method 68 , and a suitable correlation (r = 0.971) was obtained. The results were presented as μmol/gFW (Fresh weight).

Soil analyses.
The moisture content of the soil samples was determined using the weight loss method. The soil sample was weighed initially and subjected to drying at 105 ºC in an electric oven until a constant weight over time was recorded. Soil moisture content was estimated by the difference between the initial and final weights of the sample. Soil nutrient analyses were conducted only for the fine sediment component (< 128 μm). TN and TP were analyzed separately from the soil in different method. TN in the sediment (oven-dried) was analyzed by a CHN elemental analyzer (Yanako CHN coder MT-5 and Autosampler MTA-3, manufactured by Yanako CO., Ltd., Japan). TP was determined by the molybdenum blue colorimetric method 83 after digestion with H 2 SO 4 -HClO 4 (APHA, 1998) 84 .

Distribution analysis of target species.
In both river sites, all individual trees taller than 1 m were counted per species for the entire gravel bar for the Arakawa River and the area between the dikes of the Hii River along a 2.7 km long stretch 85 . Compared with the channel morphology obtained by the MLIT River Management Office of each river, the Arakawa Upstream River Management Office, and the Izumo River Management Office, tree densities were obtained for every 0.5 m elevation. Statistical analysis. Data analyses were carried out using R 86 for the dependent variables. Raw data of all variables were checked for normal distribution with the one-sample Kolmogorov-Smirnov test and for homogeneity of the variances with Levene's test. When necessary, arcsine transformation was performed. Results are presented as mean ± SD (n = 3). Data were subjected to a one-way Analysis of Variance (ANOVA) followed by Dunkan's multiple range test to evaluate the mean difference at a 0.05 significance level. The correlation coefficient (r) and the significance levels (p) were used to evaluate the strength and significance of the parameters estimated. Non-linear regression analysis was performed to view species to species effect among parameters (H 2 O 2 , soil moisture content, TN, TP, and TN: TP). A general additive model was performed to find out the efficacy of derived Monod equations among parameters (soil moisture, H 2 O 2 , TN, TP, TN: TP, and elevation).

Data availability
The authors highly appreciate and state that data is available for everyone in the supplementary file named Raw Data.